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1. Introduction 



The exploration of the phase diagram of Quantum Chromodynamics (QCD) (see, e.g., 
[1]) is a very active area of ongoing research in both, experimental and theoretical 
physics. 

It is well known that the zero temperature ground state of QCD (at normal nuclear 
densities) has two important features: chiral symmetry is broken spontaneously and the 
color charges are confined to the interior of the individual hadrons (the scale is about 
lfermi). However, the studies in lattice QCD predicted that, beyond a temperature 
of about 150MeV (in fact, the critical temperature is T c ~ 200MeV), the hadronic 
matter undergoes a transition from the confined phase to a deconfined color phase (of 
quarks and gluons). 1 

Since the transition temperature is not extremely high, it was expected that this 
phase of matter can even be produced in the laboratory. Indeed, the experimental data 
obtained at the Relativistic Heavy Ion Collider (RHIC) have confirmed the 'theoretical' 
expectation. 

What came as a big surprise, though, is that the experimental data indicate that 
the quark-gluon plasma is a new state of thermalized matter, which exhibits almost 
ideal hydrodynamic behaviour. The early thermalization and a very small viscosity to 
entropy density area (77/s) are characteristics of a strongly interacting system — at 
weak coupling (g <C 1), the equilibration time and ratio r]/s are proportional with the 
mean free path (~ 1/g 4 ) and so they are parametrically large. 

In other words, the experimental data are consistent with the interpretation of the 
new state as a strongly interacting plasma (see, e.g., [2] and references therein). This 
new phase of nuclear matter (the state of deconfinement super-hot hadronic matter) is 
known as the strongly coupled quark-gluon plasma (sQGP). 

Obviously, the perturbation theory is not useful to investigate the properties of 
sQGP. Also, the lattice QCD study is not suitable for computing the dynamical quan- 
tities of sQGP. Interestingly enough, valuable guidance for understanding the near per- 
fect fluidity observed at RHIC was obtained from the Anti-deSitter(AdS)/conformal 
field theory(CFT) duality [3]. 2 

It is important to emphasize that the gravity dual of QCD is not known. However, 
the QCD matter of interest is in a deconfined phase for which the 'conformal anomaly' 
(due to the running of the coupling) appears to be relatively small. Therefore, since 
QCD is approximately conformal at sufficiently large energies, one expects some sim- 

1 Also, by increasing the temperature and/or density the QCD vacuum restores its chiral invariance. 
2 A concise review of phenomenological problems that can not be solved by the standard field 
theoretic approach, but for which the AdS/CFT duality may be useful, can be found in [4]. 
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ilarities with some of the gauge theories that can be studied by using the AdS/CFT 
duality. 

Indeed, a study of quantum field theories (QFT) with a gravity dual revealed the 
fact that, when heated up to finite temperature, these QFTs behave hydrodynamically 
(see, e.g., [5] and references therein and, also, [6]) at large distances and time scales. 

The holographic hydrodynamics (the first attempt to study hydrodynamics via 
AdS/CFT was [7]) is an important tool for understanding some properties of strongly 
coupled quantum field theories in terms of AdS black holes physics. A notorious ex- 
ample is the computation of the ratio of the shear viscosity (if) to the entropy density 
(s). It was proposed in [8] that there is an universal viscosity bound 3 , which is also 
satisfied by sQGP: 

• (1-1) 

Kubo formula relates the shear viscosity to the two-point function of energy mo- 
mentum tensor in zero frequency limit (see Appendix A for a review). From the 
field/operator correspondence of AdS/CFT duality, we know that the energy momen- 
tum tensor of boundary field theory is sourced by bulk graviton excitations. In [7], 
the authors have considered graviton excitations polarized parallel to the black brane, 
which are moving transverse to it. They found that the graviton absorption cross 
section is related to the shear viscosity coefficient. 

In this paper, we use holographic techniques to investigate the properties of 4- 
dimensional gauge theories with finite electric charge density in the presence of a con- 
stant magnetic field. Strong magnetic fields are created in heavy ion collisions at RHIC, 
and some models were studied in [10]. 

A study of black branes solutions of Einstein-Maxwell AdS gravity with a gauge 
Chern-Simons term can be found in [11, 12]. We use the entropy function to carefully 
study the near horizon geometry of these black branes at zero temperature and discuss 
different branches with finite area horizons. This allows us to obtain an analytic ex- 
pression for the entropy density that supports the numerical analysis of [12] (though, 
we find a larger class of near horizon geometries). That is, there is a critical value of 
the magnetic field for which the entropy vanishes. 

We also extend the work of [11] by including a modulus, which is exponentially 
coupled to the gauge field kinetic term in the action. This kind of coupling appears in 
consistent truncations of type I IB supergravity on Sasaki-Einstein manifolds [15]. In 
the extremal limit, we argue that, due to the attractor mechanism, the near horizon 
geometry is stable under scalar fields fluctuations. 



3 However, there are known examples for which the bound is violated (see [9] and references therein). 
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More importantly, unlike [11], our main focus is on the hydrodynamics properties 
of this system. We use the membrane paradigm approach [16] and its higher derivative 
generalization [17] (see, also, [18]) to compute the shear viscosity (related work for 
finite chemical potential can be found in [19, 20]). In particular, we study the effect 
of the moduli (in the presence of a generic set of four-derivative interactions) on the 
viscosity bound. 

Our results are given in terms of near horizon data, and not in terms of the field 
theory parameters that would be measured at the AdS^ boundary. For a computation 
of the viscosity of our system, that is enough. However, there is a non-trivial flow 
for other transport coefficients, e.g. the conductivity. Also, it will be interesting to 
express the entropy in terms of asymptotic parameters — in [12] this problem was 
solved analytically (for vanishing scalar field). We will present some of these results in 
the forthcoming work [13, 14]. 

An overview of the paper is as follows: In Section 2, we present the bulk action, set- 
up the conventions for our model, and provide some useful holographic data. In Section 
3, we carefully study the near horizon geometry of electromagnetic black branes at zero 
and finite temperature. In the extremal limit, we use the entropy function and present 
exact near horizon geometry solutions. In Section 4, we compute the shear viscosity to 
entropy density in the presence of a generic set of four derivative interactions. Finally, 
we conclude with a discussion of our results and present some future directions. The 
appendices contain supporting material. 

2. General set-up 

The AdS/CFT duality provides a concrete relation between the regular black brane 
solutions in AdS$ and the hydrodynamic regime of strongly coupled 4-dimensional 
quantum field theories. 

In this section, we fix the conventions for our model and present a brief review of 
holography relevant to our work. We follow closely [11], though our model is different 
because we also consider scalar fields. In this way, it is straightforwardly to compare 
some of the results. However, one important difference is that we do not rescale the 
horizon radius to one, as in [11] - the physics in the presence of a scalar field is 
unambiguous if we do not use this rescaling. 

2.1 The model 

We will focus on a 5-dimensional theory of gravity coupled to a massless scalar and an 
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abelian electromagnetic field whose action is 4 

Sem = J d 5 xV=g (R+12- e a ^F, v F^ - d^cp) + S cs , 

Scs 



(2.1) 



(2.2) 

12// 2 , and also 



167TG5 

In this paper, we consider a constant moduli potential, V(<p) = 2A = — 12/Z 2 , 
fix the radius of AdS to be I — 1. 

Since the equations of motion for the gauge field simplify, we choose the constant 
in front of the Chern-Simons term to be £/3. The action (2.1) with various values for a 
resembles (truncated) actions obtained in string compactifications [15]. The coupling 
C captures the strength of the anomaly of the boundary current. 

The equations of motion for the metric, scalar, and electromagnetic field 



d u A fl ) are 



--d„ 



'-9 



9fj.u 2F^pF utJ g 
—gd^) ' 



pa 



dpipd u ip = , 



e <*eW dv {y/=jjF v ») + ^e^ 5 F pa F l5 = , 



(2.3) 
(2.4) 
(2.5) 



where we have varied the scalar and the electromagnetic field independently. The 
Bianchi identities for the gauge field are F^-x] = 0. 

Since we are interested in a theory for which the Chern-Simons term has a non- 
trivial contribution, we consider the following ansatz for the gauge field: 5 



B B 
A = E(r)dt — —ydx + —xdy — P(r)dz . 



(2.6) 



(2.7) 



The magnetic field, B, is fixed to be constant by the Bianchi identities. Thus, the field 
strength is 6 

/ E'(r) -P'(r)\ 
-E'(r) 
F= 5 
0-500 
\ P'{r) / 

4 We also have to add counterterms to regularize the action — for (2.1), the counterterms are given 
in [21] (see, also, the nice review [22] and references therein for a more detailed discussion). 

5 Please note that our P'{r) is the same as P{r) in [11] and our Z(r) in the ansatz of the metric 
(2.8) is their C(r). 

6 Our convention for the coordinates is (r, t, x, y, z). 
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where ' denotes derivatives with respect to r. 

Our analysis is on time-independent black hole solutions and so we consider the 
following ansatz for the metric 

ds * = ^f_ u(r)dt 2 + e 2V ^ (dx 2 + dy 2 ) + e 2W ^ (dz + Z{r)dtf , (2.8) 
U(r) 

which is compatible with the symmetries of the problem. 

In this case, the horizon is located at (the biggest root of) U(r h ) = and the 
temperature can be easily computed on the Euclidean section — we obtain 

T = . (2.9) 

Aix 

By using the metric ansatz (2.8), we can rewrite the Maxwell equations as 

{Q( r yV(r)+W(r)+a V (r)y _ 2(BP'(r) = 0, (2.10) 

[ e W(r)+W(r)+artr) ({/(^^Wp/^) _ Q( r ) Z (r))}' - 2(BE'{r) = (2.11) 

where 

Q(r) = E'(r) + Z(r)P'(r). (2.12) 

It is easier to work with combinations of Einstein equations rather than using 
directly (2.3). First, we extract the expressions of second derivatives of the functions 
that characterize the metric in the following way: we obtain W"(r) from (rr)- , V"(r) 
from (xx)- , and U"(r) from (^-component of Einstein equations. 

Let us now consider the (te)-component of Einstein equations in which we replace 
W"(r), V"(r), and U"(r) — we obtain 

e 2W(r) pv'(r)Z'(r) + 3W'(r)Z'(r) + Z'\r)\ - AQ{r)P\r)e a ^ r) = 0. (2.13) 

An important observation, which can be drawn by studying the system of equa- 
tions (2.10)-(2.13), is that a non-zero magnetic field is not compatible with a constant 
function Z(r) (and, also, P(r)). 

The other (independent) combinations of Einstein equations are obtained as fol- 
lows: by replacing U"(r) in the (rr)-component of Einstein equations we get 

2B 2 e a<p(r)-W(r) + 2Q( r ) 2 e «vW + U(r)[W'{r) + W {r)f 

+[U(r) (2V'(r) + W'(r))]' + X -e 2W ^Z'{r) 2 - 12 = (2.14) 
by replacing W"(r) in the (zz)-component of Einstein equations 

-AQ{rfe a ^ + U"{r) + U'{r)[2V'{r) - W'{r)] - 2e 2W{r) Z\rf 
+2U(r)[2V"(r) - 2V'(r)W'(r) + 2V'{rf + ip'(r) 2 } = (2.15) 
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and the last one is, in fact, the (xr)-component of Einstein equations 



e 2W(r)[_ 4S 2 e a V (r) _ Ze^U' {r)V {r) + 12e 4y(r) ] - 2Q(r) 2 e AV{r)+2W{r)+aip{r) 
+U(r)e 4V ^[2P'(r) 2 e a ^ - 3e 2W(r ^ (V"(r) + V'(r)W(r) + 2V'{r) 2 )] = 0. (2.16) 

We also use the ansatz of the metric in the equation of motion for the scalar (2.4) and 
so this equation becomes 

[U( r )e 2V W +w W<S(r)Y + ae 2V ^ +w ^ +a ^ [B 2 e~ 4V ^ + U{r)e~ 2W ^ P'{rf - Q{r) 2 ] = 0. 

(2.17) 

Due to the non-trivial coupling between the scalar and gauge fields, the equation 
(2.4) has a non-trivial right hand side. The non-vanishing electromagnetic field may 
also be understood as a source for the scalar field. Thus, the scalar charge is determined 
by the electric and magnetic charges and so it is not an independent parameter that 
characterizes the system — this charge plays an important role when the asymptotic 
value of the scalar is not fixed (see [24]). 

We would like to conclude this section with an observation on the Hamiltonian 
constraint. A vanishing Hamiltonian is a characteristic feature of any theory that is 
invariant under arbitrary coordinate transformations — for our system, we can obtain 
a first order differential equation by replacing W"(r),V"(r), and U"(r) in the (tt)- 
component of Einstein equations. 

The Hamiltonian constraint, which can be enforced as an initial condition, has the 
following expression: 

2B 2 e a<p(r)-4V(r) + jj ^ [_ 2 p'( r )2 e a V (r)-2W(r) + ^ V '(r)W'(r) + 2V» 2 - (2.18) 

-V?'(r) 2 ] + 2Q(r) 2 e a ^ + 2U'(r)V'(r) + U'(r)W'(r) + - e 2W ^Z'(r) 2 - 12 = . 

2 

2.2 Holographic data 

The AdS/CFT correspondence [3] is a concrete example of open/closed string duality. 
Remarkably, at strong coupling M = 4 SYM is best described by classical supergravity 
— in this section, we closely follow the reviews [22, 25]. 

The observables in the field theory side of the duality are the correlation functions 
of gauge invariant operators, which are composites of the elementary fields. Any su- 
pergravity field, $, corresponds to an operator, O, in the (boundary) field theory — in 
particular, bulk gauge fields correspond to boundary symmetry currents. 

We are interested in planar black hole solutions in AdS and so we will consider the 
topology of the AdS timelike boundary to be R 3 x R + — the intrinsic coordinates of the 
boundary are denoted by x and the radial coordinate by r. A concrete prescription to 
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do computations on the Euclidean section was given in [26]. That is, the on-shell bulk 
partition function with the boundary conditions $(x, r)^ = $o(x) is the generating 
functional of the boundary QFT correlation functions: 

Zsuqra[*o{x)\ = [ / e J***o(*>°(*>\ . (2.19) 

J \ I QFT 

We work in the saddle point approximation and so the generating function of QFT is 
determined by the on-shell supergravity action, Wqft[^o] — —S nsheii[&o\- To compute 
the two-point function of a particular operator in the field theory, we should first solve 
the linearized equation of the corresponding supergravity field with the appropriate 
boundary conditions, then evaluate the on-shell sugra action, and, finally, compute the 
second functional derivative with respect to the source: 

< 0(x)0(y) >= - J) S ZTf S • ( 2 - 2 °) 

We are interested in a bulk gravity theory with an electromagnetic field non- 
minimally coupled to a scalar. In general there are two types of perturbations of 
AdS: those that modify the bulk but preserve the AdS asymptotics (e.g., a black hole 
is interpreted as a finite temperature state in dual QFT) and those that do not modify 
the bulk and act as sources for operators in the dual QFT. The AdS metric (graviton) 
couples to the stress tensor and the gauge field couples to the R-charged current 
of the field theory. 

To obtain the viscosity (conductivity) we should compute the two-point functions 
of the stress tensor (current). Due to the coupling of the gauge field with a scalar, there 
will be a non-trivial 3-point function mixing two currents with a scalar operator — in 
this work, we are not interested in this correlation function. 

The bulk solution is a geometry with a horizon (at zero or finite temperature) 
and to get the information about the hydrodynamics one should compute the retarded 
Green functions directly on the Lorentzian section. One reason is that, on the Euclidean 
section, there is a discrete set of frequencies, the Matsubara frequencies (u = 2irinT 
for bosons and u = mT{2n + l)/2 for fermions). However, since we are interested 
in a small u in the hydrodynamic approximation, an analytic continuation from the 
Euclidean section where there is a minimum frequency (u m = 2niT) is not useful. 

The general prescription for real-time retarded Green's functions within AdS/CFT 
duality can be found in [27, 28]. The existence of hydrodynamics modes is reflected 
by the existence of the poles of the retarded correlators. However, to compute the 
viscosity bound, we will use a different prescription, namely the membrane paradigm 
proposal of [16]. 
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We would like to present now a discussion on the boundary conditions for gauge 
fields. The bulk gauge field (2.7) is invariant under translations in the boundary coor- 
dinates (x). The component tangent to the boundary of the gauge field is constant at 
large r, which corresponds to a magnetic field B z perpendicular on the (xy)-plane. 

If there is an electric field, E x , in x-direction, then the effect of a magnetic field in z 
direction induces a current along ^-direction. For, A t (x) 7^ and Ai = A y (x) 7^ 0, the 
conserved currents on the boundary are the charge density, q, and the current density, 
J — they are defined as 



5S 
q= 6A t 



r 63 



(2.21) 

boundary 



boundary SA. 

To obtain the boundary current, we have to impose the following boundary conditions: 

3Af I boundar y = , 5 A j | boundary free. (2.22) 

The condition of regularity of the gauge potential at the horizon requires that A t 
should vanish (on the Euclidean section, the radius of the Euclidean time shrinks to zero 
at the horizon). Therefore, we should add a gauge term, T so that A t {r h ) — T = 0. This 
term, which remains constant at the boundary, plays the role of a chemical potential 
for the electric charge. 

Let us end up this section with a comment about the Chern-Simons term. The 
coupling of the gauge Chern-Simons term is proportional to the chiral anomaly of the 
corresponding current in the dual field theory. The current anomaly is given by a 
one-loop triangle Feynman diagram to the 3-point function of (R-) currents. Therefore, 
there also exists a chiral magnetic current which is induced parallel to the applied 
magnetic field, J z ~ (BY. 



3. Near horizon geometry 

Shear viscosity of the boundary fluid can be computed from the knowledge of the near 
horizon physics only [16] - in Section 4, we will present a detailed analysis. This 
observation is very useful when the bulk Lagrangian is very complicated. For example, 
for the model we are interested in, gravity is coupled to various matter fields in a 
non-trivial way and it is technically difficult to find an analytic solution of the system. 

However, it is not a very complicated problem to find the near horizon geometry, 
i.e. how the metric and other fields behave in the near horizon limit. The reason is 
that, in principle, we do not have to solve any differential equation to find the near 
horizon geometry. A study of the extremal near horizon geometries is interesting in 
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its own because it provides information about the instabilities that may appear in the 
theory. 

The algorithm is as follows: first we obtain the field equations. Then, we consider 
a suitable near horizon ansatz for different fields — all fields should be regular at the 
horizon. Substituting the ansatze for fields in the corresponding equations of motion, 
one can solve them consistently order by order in (r — 77J. 

In this section, we will first find the near-horizon geometry of the non-extremal 
black hole, which we will need to compute the shear viscosity to entropy density ratio 
in Section 4. In the extremal limit, we will see that, due to the attractor mechanism, 
the near horizon geometry is universal regardless of the asymptotic values of the scalars. 
We present a detailed analysis of the branches of solutions with finite horizon area. 

3.1 Non-extremal case 

As in [11], we work with a coordinate system in which the solution takes the canonical 
form at the horizon. That is, the field strength F H and the metric ds 2 H are 



where q and B are the charge density (of the black brane) and the magnetic field at 
the horizon, respectively. In this way, the gauge freedom is removed and the initial 
conditions are 



A similar analysis (and numerical solutions) in the presence of the Gauss-Bonnet term 
but without the Chern-Simons term was presented in [29]. 

The generic solutions have a non-degenerate horizon. Near the event horizon, they 
admit a power series expansion of the form (using the definition of the temperature 
(2.9) in the expression of U): 



U (r) = AnT{r - r h ) + u 2 (r - r h ) 2 H , 

V(r) = ln(r h ) +v 1 (r- r h ) + v 2 (r - r h ) 2 H , 

W{r) = ln(r h ) +w ± (r- r h ) + w 2 (r - r h ) 2 + ■■■ , 

Z{r) = zi(r - r h ) + z 2 {r - r h f H , 

E{r) = q(r - r h ) + qi (r - r h f + q 2 (r - r h f + • • • , 
P(r) = p( r - r h ) +p 1 {r- r h ) 2 H , 

(p(r) = ip h + y?i(r - r h ) + ip 2 {r - r h ) 2 H (3.3) 



Fh = q dr A dt + B dx A dy — p dr A dz, 
ds 2 H = r 2 H {dx 2 + dy 2 + dz 2 ), 



(3.1) 



U{r h ) = Z(r h ) = P{r h ) = , V(r h ) = W(r h ) = \n(r h ). 



(3.2) 
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It is important to emphasize that, what is generally called near horizon geometry for a 
non-extremal black hole is just a truncation of the above series expansion. To compute 
the shear viscosity, though, we need also some data at the order (r — r h ) 2 . 

Another observation is that, in principle, one can use a boost transformation in z 
direction to set p = (see [11]). However, the boost transformation is singular at some 
point outside the black hole horizon. In our analysis we keep the value of p non-zero 
and determine it in terms of other horizon data. We will see in Section 4 that the 
expressions for the entropy, shear viscosity, and their ratio remain unchanged if we set 
the horizon value of P'(r) to be zero (in other words, they do not depend of p). This 
is expected due to the fact that the physical quantities should be invariant under the 
boost transformations. 

By substituting the ansatz (3.3) in the field equations, we get the following expres- 
sions for the coefficients at the order (r — r h ): 7 

2B 2 e aVh + r{ (q 2 e°"? h - 6) 



Wi 



AB 2 e avn _ r 4 ^ q 2 e a<p h + 3^2 _ 34) 



p 



q (2B(e- a ^ + z x rl) 
4nTr h 

-2aip h / 

^ = ^l(25 2 [(« 2 + 2)e— +4C 2 ] 

_ r 4 e 2a W ( a 2 _ 2 ) e a Vh + ^ + 2 4]^ . (3.4) 

However, in higher derivative gravity theories the only coefficient at the order 
(r — rh) 2 , which we need for viscosity bound computation, is u 2 . But, for completeness, 
we present the expressions of all the other coefficients that appear at order (r — r h ) 2 , 
in Appendix B. 

A non-extremal charged scalar black hole is characterized by four independent 
parameters: the mass, electric charge, magnetic field, and also the value of the scalar 
at the horizon, (p^. In this case, the horizon radius (and so the entropy) and the horizon 
value of the scalar depend of the asymptotic boundary data (ty?oo)- We will see in the 
next subsection that this is in contrast with the extremal case for which we obtain an 
attractor behaviour at the horizon. 



7 We obtain the results as functions of the coefficients (T, q, B 7 ifh, rh, and z\) — this will simplify 
the computations of the shear viscosity. 
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At first sight, it may seem surprising that the data (3.4, B.l) we need to compute 
the entropy and shear viscosity depend also on zi, a coefficient that we do not com- 
pute explicitly. However, we will see in Section 4 that the final values of the physical 
quantities depend in fact just on four independent parameters, namely (q, B,rh,<fh) 
- we 'trade' the mass for the horizon radius and so the independent parameters that 
completely characterize the black hole are the ones mentioned above. 

3.2 Extremal limit and attractor mechanism 

It is by now well known that the extremal black holes in theories of Einstein gravity with 
scalar fields non-trivially coupled with abelian gauge fields have an enhanced symmetry 
of the near horizon geometry [30]. That is, the near horizon geometry contains an AdS% 
spacetime. 

The attractor mechanism [31] is also valid for non-BPS extremal black holes [32, 
33, 34] 8 — a qualitative explanation is that, due to the infinite long throat of AdS%, 
the near horizon geometry has no memory of asymptotic data. Consequently, the near 
horizon geometry is universal and the entropy does not depend on the asymptotic values 
of the scalars. 

In this section we present a careful analysis of the extremal near horizon geometry 
by using directly the equations of motion and also the entropy function formalism. The 
advantage of the latter is that we can compute the physical charges and obtain analytic 
expressions for the entropy density. In this way, we can confirm the numerical analysis 
of [12]. 

The (non-supersymmetric) attractor mechanism in AdS and its embedding in string 
theory were studied in [37]. In the rest of this section we follow closely [37, 29], though 
we should keep in mind that our goal is not to investigate just the flow of the entropy 
but also the flow of some hydrodynamic transport coefficients, and we also consider a 
non-trivial Chern-Simons term. 

3.2.1 Attractor mechanism 

In AdS spacetime, the BPS condition is different than the extremal limit. For static 
black holes, the BPS limit is a naked singularity, but in the extremal limit the entropy 
can remain finite. 

Unlike the non-extremal case, in the extremal case the near horizon geometry is an 
exact solution of the equations of motion and not a truncation in a Taylor expansion. 

8 A generalization to stationary non-susy black holes was given in [35]. More details on the non- 
supersymmetric attractors can be found in the nice reviews [36]. 
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Let us consider the most general ansatz for the near horizon geometry: 

/ dr 2 \ 

ds 2 = L ( — - r 2 dt 2 J + vi (dx 2 + dy 2 ) + v 2 (dz + z x rd£) 2 , 

F = qdr A dt + Bdx A dy, 
ip(r) = ip h . (3.5) 

Another important difference with the non-extremal case is that, due to the enhanced 
symmetry of the near horizon geometry, p vanishes (the z component of the gauge field 
is constant). As a consistency check, it can be shown that with this ansatz, in the 
absence of the scalar field, we obtain the results of [ll]. 9 

Let us consider now in more detail the non-trivial case when the scalar is turned 
on. The first Maxwell equation can be easily integrated and we obtain 

_ v 2 v 2 q 2 e 2a ^ 

where the integration constant is related to the physical charge. 10 
We use the following system of independent equations: 11 

4B 2 L 2 v 2 e 2a ^ h + v\v 2 e^ h (24L 2 - 6L + Zv 2 z() + 8L 2 Q 2 p = 0, 

4B 2 L 2 v 2 e 2a ^ h - 4L 2 (Q 2 p - Qv\v 2 e aiph ) - ?,v{v 2 2 z\e a ^ = 0, 

Q2 -acp h 

_2 B * e <*PH _ ^II + = 0) 

v 2 

V2 ^ +24 y^tp_^^ (3 . 7) 

2B 2 v 2 e a ^ _ v x v 2 2 z\ 

Qp[ + 4zi =0, 

a [ B 2 L 2 - ) = 0. 

v 2 



9 We emphasize again that the authors of [11] have shifted the radial coordinate, r, to set the horizon 
at r h = 1. 

10 Q P includes the contribution from constant z component of gauge field. 

11 The first one is the (rr )-component of Einstein equation, the second equation is a combination of 
(rr)- and (tt)-components, the third one is the (xx)-component, the fourth one is the (te)-component, 
the next one is a combination of (xx)- and (zz)-components, and the last two are the remaining 
non-trivial Maxwell equation and the scalar equation, respectively. 
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As we already explained above, we expect that the near horizon geometry to be 
completely fixed by the charges. By solving the system of equations (3.7), we obtain 
the following expressions for AdS 2 radius, size of transverse space, strength of fibration, 
and horizon value of the scalar: 



work with planar black holes — for extremal spherical black holes in AdSs, the radius 
of AdS 2 does not match the radius of AdS 5 [37]. 

One important observation is that the horizon value of the scalar is fixed by the 
Chern-Simons coupling. We will reobtain this result by using the entropy function. 

For a regular extremal solution whose entropy does not vanish, there is a large 
ground state degeneracy. However, an important question is then, if the entropy is 
'stable' against changes of the boundary conditions for fields in the bulk. In [11], the 
authors argue that the near horizon geometry (AdS 2 x R 3 ) of the purely electrically 
charged brane is unstable under the addition of a magnetic field. 

In the rest of this section we argue that, in the presence of the scalar and electro- 
magnetic field, the entropy does not change under general perturbations of the scalar. 
A similar discussion for a 5-dimensional theory without the Chern-Simons term can be 
found in [37, 29]. 12 

The arguments are as follows. Due to the non-trivial coupling between the scalar 
and gauge fields, there is an 'effective potential' for the scalar [34]. If the effective 
potential has a stable minimum at the horizon, the solution is regular and the entropy 
is determined completely by the charge parameters (if there are flat directions, not all 
moduli are stabilized, but the entropy does not depend of the flat directions). 

We can see that, without a gauge field, the theory is scaling invariant and so the 
scalar field is always massless. By turning on just one gauge field (electric or magnetic) 
the symmetry is broken and one expects that the solution is not regular. The reason is 
that, even if an effective potential is generated, the potential has no stable minimum. 

However, by turning on a second gauge field (in our case, we consider an electro- 
magnetic field) we expect to obtain regular solutions. The effective potential can have 
stable minima and the near horizon geometry is universal. 

One important observation, though, is that the ansatz of our solution is different 
than the one of [34] - we have non-trivial terms (dzdt) in the metric. A concrete 
expression for the effective potential in this case is not known. However, in the near 

12 Recently, the authors of [38] have used the same kind of arguments to discuss charged scalar black 
holes in AdS±. 




The radius of AdS 2 is the same as the radius of AdS$. This is due to the fact that we 
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horizon limit the effective potential method should be equivalent with the entropy 
function method of Sen. In what follows, we will use the attractor equations to find 
the near horizon geometry (and so the entropy density) of the extremal solutions. 

3.2.2 Entropy function 

By using the entropy function formalism we can explicitly compute the physical charges 
and obtain an analytic expression for the entropy density. Since our theory contains 
the gauge Chern-Simons term, which is not gauge invariant, we have to Kaluza-Klein 
(KK) reduce our action to four dimensions. 

Let us consider first the case with the scalar turned off. We consider the z direction 
compactified on a circle of radius j3 and a general KK ansatz 

g a pdx a dx p = G ab dx a dx b + G AB {dy A + A A dx a )[dy B + A B dx a ), 

A (5) = A (5) dx n = Afdx a + C B {x a ){dy B + A B dx a ), (3.9) 

where a, b are 4D indices and A, B are compact indices [z in our case). With this nota- 
tion we have splitted the coordinates as x M = (x a ,y A ) and so Affi is the 5-dimensional 
gauge potential, A^ is the 4-dimensional gauge potential, and A B is the KK gauge 
potential. 

We use the following results of the dimensional reduction [23] (there is no depen- 
dence of KK coordinates): 

= v /Z G V / det(G AB ) , F^ 2 = F^F® = F^F® + 2F^ aB F^ , 

d D s~iacs~<bd/~< ipAfpB , ^ o s~< f\a/~iAB ^/iiBa fi s~iCDp\as~< 

K5 — -K4 — -Lr Lr LrAB-fab-^cd + ^° a< ~' AB ° ^ ~~ 4^ Ca^AB^ O LrQD ~ 

-d a {G AB d a G AB ). (3.10) 

Now, we are ready to write down the ansatz we are interested in — the metric, gauge 
field, and relations between the 5-dimensional gauge potential, KK gauge potential, 
and 4-dimensional gauge potential are 

(if 



ds 2 = L{-— + r 2 dr 2 ) + vi{dx 2 + dy 2 ) + v 2 (dz + z^rdt) 2 , 



B B 

A^dx* 1 = firdt — —ydx + —xdy — p\{dz + z±rdt), 



2 a 2 

Fif =9 = 0- PlZl , = F« = B , Flf = , F* t = z x . (3.11) 
The on-shell action and entropy function are 

A X yP 



s 



IQnG 



5 



f-^-f 2B 2 2{d- PlZl ) 2 v 2 z 2 -4L \ , PlZl \ 

{-^f + V + —21?- + 12 ) + WPlC v~ " V 



£ = 2nPA xy [Qd + Q Zl - — — ) , (3.12) 

P-^-xy , 
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where Q is the 4- dimensional physical charge, is the physical charge associated to 
KK gauge field, and A xy = f dxdy. 

The equations of motion in the near horizon limit are 



77 x/i 2 [AB 2 L 2 + v 2 (-AL + 24L 2 + Aq 2 + v 2 zf)] = 0, 



Lv( 
3ti 



[AB 2 L 2 - vj(-AL + 2AL 2 + Aq 2 + 3v 2 zj)] = 0, 

[AB 2 L 2 + vj(-2AL 2 + Aq 2 + v 2 zj)] = (3.13) 

^ [L (O - ABpjC) + ^^(^ig - v 2 zi)] = 0, 

^ [2BL( + v ly fU- 2Zl ] = 0, 

2vr/3 ^8BpiC - 4t;i ^ 9 + <2j = 0. 

First, we solve the last equation to obtain the 5-dimensional electric potential. Then, 
we use this result in all the other equations, and from the forth equation we obtain the 
following expression for the KK potential: 13 

L( Pl Q + 9 + ABpjQ 
Zl = y~2 • ( 3 - 14 ) 

ViV 2 ' 

Next, we compute pi from the sum of the first two equations — we obtain a 4th order 
equation with the following four solutions: 

Q + y/Q 2 - 16B(2Bv 2 + 9K -Q + ^Q 2 - l6B(2Bv 2 + gg 
P 8B( ' P SB( 

-Q + VQ 2 + l6B(2Bv 2 - 9"K Q + VQ 2 + l6B(2Bv 2 - 9~K 

P 8B( ' P 8B( 

In the absence of the magnetic field, we expect a trivial fibration for our geometry: we 
can see that the physical solutions are the second and third ones because p\ is finite in 
the limit B — > 0. Therefore, depending on p\ solution, we have two distinct families of 
solutions. To compute the other near horizon parameters, we proceed as follows: by 
subtracting the third equation from the first one we solve for v 2 , and get 

V2 ~ 8[v 2 + AB*L(-l±()Y (3 - 15) 



13 



We note that by combining the first equations we get L = 1/12. 
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The +/— solution corresponds to the second and third solution for p 1 , respectively. 

The 5th equation, which is very important for understanding the solutions, and 
the entropy are 



It is now clear that, from the equation (3.16), there are three distinct branches: 14 
B = 0, q = 0(v 2 l = B 2 /3), ( = ±1. (3.17) 
From (3.16), we can also read off the corresponding entropy for all three branches: 

Branch 1: B = 

This branch corresponds to the usual Reissner-Nordstrom solution and its entropy 

is 

S = P A *y*Q . (3.18) 

Branch 2: q = [Bjv x = yg) 

We compute all near horizon parameters in terms of physical charges. We obtain 
Pi = — g§^ and also find two possible solutions for z\ and v 2 : 



^V2B 1 / Q 2 \ 

setl : zi = , v 2 = — — — - 16B0 , (3. 

_ 16B0 32B V C / 



19) 



or 



set2 : zi = 
The corresponding entropies are 



2^6 1 / Q 2 \ 

~ / > ^ = 7^ 16B0 - — . (3.20) 

x/3y 16B6 - ^ 32B V C / 



5l = J^! _ i 6se , , S 2 = £^Ji6£6-^. (3.21) 



14 We rewrite g in terms of the physical parameters g = V ^^^' B - to obtain the constraint between 
V\ and B for g = 0. 
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Branch 3: ( = ±1 



In this case, one of the near horizon equations of motion, in fact, is used to get a 
fixed value for (. Therefore, we can not compute all the parameters, one of them can 
not be fixed. However, as expected, the entropy depends only on the physical charges: 



Now, we are ready to compare our analytic results with the numerical results of 
[12]. We just saw that with the entropy function formalism we get a similar condition 
as the one of [11, 12], namely qB(l ± Q — 0. Therefore, we also obtain solutions for 
qB — or ( — ±1. In the first case, when ( ^ ±1, there is a regular solution (when 
B = 0, branch 1) with the expected near horizon geometry (AdS2 x i? 3 ). 

For the case q = (branch 2), depending on the sign of (, we have different 
situations. Let us assume that the physical KK charge is O > 0. When ( > (but 
7^ 1), the positivity of t; 2 (which is also reality of Z\) implies that B has a critical 
value B c = y^L_ (solution set 1). Once B crosses that critical value, the solution set 
1 becomes singular and the solution set 2 takes over. In this case, for a generic value 
of the magnetic field, the near-horizon geometry is always AdS 3 x R 2 (see, also, [11]). 
For C < (and ^ — 1), the solution set 1 does not exist for any positive value of the 
magnetic field. In this case, the only possible solution is set 2, which is a regular (finite 
horizon area) solution for any positive value of B. We would also like to point out that, 
when the Chern-Simons coupling vanishes, it seems there is no finite area solution in 
the case q = and non-zero magnetic field. 

In the second case, ( = ±1, there is a family of solutions for which 2{B/vi) 2 + 
(q/L) 2 — 6 = (this constraint is obtained by solving the third and fifth equation in 
(3.13) and matches the condition of [11] where the horizon radius is rescaled so that 
Vi = 1). In this case for ( — 1 there exists a critical magnetic field B c = ^ above 
which the solution becomes a naked singularity. For ( = —1 the solution exists for any 
value of B. 

The near horizon geometry smoothly interpolates between AdS2 x R 3 (for B = 0) 
and AdSs x R 2 (for q = 0) — in between, the near horizon geometry is in fact an 
warped AdS 3 x R 2 . 15 

Let us know discuss the case when the modulus is turned on. The analysis is very 
similar with our computations above — we present the details in Appendix C. Unlike in 
the previous section where we have solved the equations of motion in the near horizon 

15 The AdSz spacetime is a special fibration of S 1 over AdS2 and the case q — corresponds to this 
special fibration. 




when C, = ±1. 



(3.22) 
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geometry, by using the entropy function we can compute the physical charges in four 
dimensions. In five dimensions, due to the existence of the gauge Chern-Simons term, 
one has to be careful about the definition of the physical charge — we will present a 
detailed analysis of this aspect in the forthcoming work [13]. 

By using the attractor equations from Appendix C, we obtain the same horizon 
value of the modulus 

Vh = — In |C| (3.23) 

as in the previous subsection. In general, due to the attractor mechanism, in the 
extremal case the horizon values of the moduli are fixed by the charges. Interestingly 
enough, this is not the case here: the modulus is fixed by the Chern-Simons coupling 
alone. We will present in Discussion section an argument of why this is the case. 

Another observation is that the physical magnetic field is also fixed by the Chern- 
Simons coupling: 

B 

— = ±v / 2e" W2 (3.24) 

Vi 

or, equivalently, the parameter v\ is fixed only by the magnetic field parameter and the 
Chern-Simons coupling and does not depend on the electric field. This is due to the 
fact that we can rescale x and y coordinates in our solution to make v\ — 1, and this 
rescaling affects just the magnetic field — that can be easily seen from our ansatz for 
the metric and gauge potential (3.11). 

In Appendix C, we present a detailed computation of the entropy — the result is 



Q 2 - IQBQ( r , 9( .x 
S = PA xy nJ . (3.25) 

It is not our goal to make a detailed analysis with the scalars turned on here (for 
different a and 'anisotropic near horizon geometries'), but we observe that, when just 
the electric or magnetic field is turned on, the solution is a naked singularity. This 
is similar with our discussion at the end of the previous subsection, when without 
the Chern-Simons term, we have argued that the effective potential does not have a 
minimum and so there is no regular solution [37, 29]. We do not know the form of the 
effective potential when the Chern-Simons term is present, but that can be easily seen 
from the equation of motion of the scalar (the last equation in (C.l) ). Consequently, 
there is just one branch of finite area solutions when the modulus is turned on. 

As in the case with the scalar turned off, there also exists a critical value of the 
magnetic field for which the entropy shrinks to zero (for ( > 0). So, the geometry 
is stable against scalar fluctuations but not against the magnetic field. However, this 
interpretation (and also the interpretations of [11, 12]) should be taken with caution. 
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It is expected that the higher derivative corrections will 'dress' the singularity with a 
horizon and so a more general analysis is important in this context. 

4. Shear Viscosity to Entropy Density Ratio 

For black objects with translation invariant horizon, for example black brane geometry, 
one can also discuss the hydrodynamics — long wave length deviation (low frequency 
fluctuation) from thermal equilibrium. 

In addition to the thermodynamic quantities the black brane is also characterized 
by the hydrodynamic parameters like viscosity, diffusion constant, electrical conduc- 
tivity, et cetera. The black _D3-brane geometry with low energy fluctuations (i.e. with 
hydrodynamic behaviour) is dual to a finite temperature gauge theory plasma living 
on boundary with hydrodynamic fluctuations. 

In this section we will study the low frequency behaviour of boundary plasma dual 
to the gravity model presented in Section 2. To do this, we use the membrane paradigm 
proposal of [16] and its generalization to higher derivative gravity theories [17]. 

4.1 Membrane paradigm and viscosity bound 

To an external observer, the black hole physics appears to be equivalent to the physics 
of a dynamical fluid membrane [39, 40]. In other words, the black hole is equivalent 
with a set of surface charges and currents at the stretched horizon. 16 

In the Lorentzian prescription of [27] for computing the Green functions, one should 
impose in-falling boundary conditions for the fields at the horizon. However, when 
computing the two-point functions, just the contribution from the boundary has to 
be taken. An explanation of why the surface terms coming from the horizon must be 
dropped was given in [28] (see, also, [41]): the boundary conditions should be defined 
in Kruskal coordinates and one has to work with the full analytic extension of the black 
hole. 

An important observation, which relates the work of [27] with the membrane 
paradigm, was made in [16] (see, also, [42, 18]). That is, the response functions in 
AdS/CFT are similar with the membrane response when the membrane is 'pushed' 

16 The stretched horizon is a timclike surface just outside the true horizon. We emphasize that the 
charges and currents are fictitious in the sense that a falling observer through the stretched horizon 
does not detect any surface sources. However, for an external observer, their presence is consistent 
with all external field configurations. 
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at the boundary. In this way, the membrane paradigm physical quantities are in fact 
concrete gauge theory observables. 17 

Interestingly, it was observed in [43] that only with the knowledge of the near 
horizon geometry one can easily calculate the shear viscosity of boundary fluid. One 
does not need to know the full analytic solutions of Einstein equations — this method 
is especially useful in higher derivative AdS gravity theories. 

As a warm-up exercise, let us consider first the model presented in Section 2.1. We 
apply the method of [16] to compute the shear viscosity coefficient of the boundary 
fluid for the action (2.1). 

Let us consider a metric perturbation of the form: 

9*v = 9 { S + M r ' x ) = & + eHr,x)}. (4.1) 

At first sight, it seems that the dual gravitational mode (4.1) does not generally 
decouple. Interestingly enough, the decoupling occurs when the momentum vanishes 
and this is what we need for the computation of the viscosity in the hydrodynamics 
limit — we provide a detailed derivation of this claim in Appendix D. 18 

By plugging (4.1) in the action and keeping the terms at order e 2 (at the first order 
in e, we obtain the equations of motion for gravitons), we get the following effective 
action for the perturbation: 19 

S = y^r / j0pdr J2 AAr, k)<t><*\r, -k)^\r, k). (4.2) 

Here, we use the Fourier transform to work in the momentum space k = {—u, k} 

*(r, x) = J 704^-0(r, k) , 0(r, —k) = 0*(r, k) (4.3) 

and 0( p )(r, k) denotes the p th derivative of the field <f)(r, k) with respect to r (p + q < 2). 

Next, we integrate by parts to obtain the bulk action for the graviton in the fol- 
lowing form (up to some total derivative terms): 

S = TTrV / 4\z dr lMr, k)<f>'(r, k)<j>'(r, —k) + A (r, k)cf>(r, k)J>(r, -k)], (4.4) 
lb7rG 5 J {2np 

17 The infalling boundary conditions at the horizons are the regularity conditions in the membrane 
paradigm: near the horizon the fields depend trough a non-singular combination of r and t, namely 
dv = dt+ ^/grr/gttdr. 

18 We would like to thank Alex Buchel for a discussion on this point. 

19 The terms that contain the derivatives with respect to the spatial coordinates, x, combine in terms 
whose coefficient is proportional with p 2, . Since we work in the hydrodynamic approximation p = 0, 
these terms do not play any role in our analysis. 
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where 

1 p 2V(r)+W(r),,2 

Mr, k) = --e^ +w ^U{r), Mr, k) = - ■ (4.5) 

At this point, it is important to emphasize that there are many total derivatives in this 
action that do not affect the equations of motion for the graviton. For the computation 
of the imaginary part of the two-point function, the coefficient of the term in 
the bulk action is important. The other total derivatives in the bulk action and the 
Gibbons- Hawking surface term contribution exactly cancel on the boundary [18]. It 
was also shown [18] (it is straightforwardly to check it also in our case) that, in the 
case of Einstein gravity, the ratio of viscosity and entropy density is not affected when 
the matter fields are minimally coupled. The effective coupling [45, 18] is 

1 A 1 (r,k) 1 

and so the viscosity coefficient of the boundary fluid stress tensor is 

V = e w ™ +w ™(-2K eff (r h )). (4.7) 

In this case, the shear viscosity to entropy density ratio turns out to be universal, 
namely 

77 1 



(4-8) 

S 4:71 



4.2 Four derivative action 



Let us now consider the action (2.1) supplemented with the most general four-derivative 
interactions [19]: 



a' 

Sud = S E m + — „ / d b xy/^g 
167rG ; " 



Cl R abcd R abcd + c 2 R abcd F ab F cd + c 3 (F 2 ) 2 (4.9) 
+c 4 F 4 + c 5 e abcde A a R bcfg R d J 9 



Since in supergravity actions the gauge kinetic terms couple to various scalars, it 
will be interesting to understand the role of the moduli in computing the viscosity 
bound. Unlike [19], our action contains a scalar, </?, and the coefficients q depend on 
the value of tp. This resembles the four-derivative supergravity action [46]. 

We treat the higher derivative terms perturbatively and apply the method of [17] 
to compute the shear viscosity coefficient of the boundary fluid. However, to obtain 
the viscosity bound we also need the entropy density. We start by using the Noether 
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charge formalism of Wald [47] (see, also, [48, 29] for a discussion in AdS) to compute 
the entropy density — we will need just the data in Section 3.1 and Appendix B. 

When we add higher derivative corrections to the action, the entropy is no longer 
given by the area law — instead, we use a general formula proposed by Wald 

f dL 

s = -2tt / — e ab e cd , (4.10) 

where e a b is the binormal to the surface %. 

By using (4.10), we obtain the following expression for the entropy density: 

s = + ^ r 3 h [d(3*M - 4«2) - 2c 2 g 2 ] + 0(a' 2 ). (4.11) 
We use the expression of u 2 given in Appendix B to rewrite this expression as 



s = 



3 „,/ ^,3 

12 



4G, G H 



3 V r h / 2 



+ 0(a' 2 ). (4.12) 



As expected, the entropy density depends on four independent parameters, namely 
(rh,q, B,tph)- Since in Wald formula only the four derivative interactions that involve 
the curvature tensor are important, the entropy only depends on c\ and c 2 (c 5 does not 
appear because the binormal has just rt components and the contribution from this 
term vanishes). 

To compute the four derivative corrections to the shear viscosity coefficient, we 
have to find the quadratic action for the transverse graviton moving in the background 
spacetime. As in previous section, we consider again the following metric perturbation 

9xy = 9 { $ + h xy (r,x) = gWj[l + e$(r,x)], (4.13) 

where e is an order counting parameter. 

In the presence of a generic n derivative term in the bulk Lagrangian, the action 
(in momentum space) can be written as 

S = / —dr A P , q (r, k)<P ip \r, -k)^\r, k). (4.14) 

5 p,q=0 

Here, 0^(r, k) denotes the p th derivative of the field 0(r, k) with respect to r and 
p + q < n. The coefficients A Ptq (r, k) depend, in general, on the coupling constant a'. 

It is important to emphasize that A Ptq with p + q > 3 are proportional to a' and 
vanish in a' — > limit [17] — these terms appear as an effect of higher derivative terms 
in the action. 
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This action does not have the canonical form as in the two derivative case. There- 
fore, it is not obvious how to generalize this approach for higher derivative case. The 
proof given in [16] was based on the canonical form of gravitons action. This problem 
was solved in [17] and we use this method to compute the shear viscosity coefficient. 

Let us now write the effective action for the transverse graviton in the canonical 
form with arbitrary coefficients: 



16nG 



1 r dud 



dud 3 k 
7 



dr 



Si(r, k)<f>'{r, -k)<jf(r, k) + B (r, k)<f>{r, k)<f>{r, -k) 



(4.15) 



We demand that the equations of motion obtained from the action (4.15) and 
original action (4.14) match at order a'. By comparing the equations of motion for 
0(r, k) from the two actions, we get the function £>o and B\. 

Once we have the effective action for 0(r, k) in the canonical form, the effective 
coupling can be easily read off: 



K ef f(r) 



1 Bijr, k) 
167rG 5 J=gg rr 



(4.16) 



where g rr is the rr-component of the inverse perturbed metric and yj—g is the deter- 
minant of the perturbed metric. Therefore, the shear viscosity coefficient is 



V = ri[-2K eff (r = r h )). 



(4.17) 



Evaluating the effective coupling in the near horizon we obtain the shear viscosity 
coefficient 



V = 



16tt G 5 

, a' ( Cl r 4 h (8q 2 e a ^ 



1' h 



327TTW! - 36vrTwi - 10u 2 + 48) - B 2 (c 2 - 8c ie a% )) 



87iG 5 r A h 



+0(a' 2 ) 

which can be rewritten as (we use the results in Appendix B) 



a r 



V = 



167rG B 2ttG, 



ci q 



o ai Ph 



+ 



4 < 



+ 0(a> 2 ). 



The ratio of the shear viscosity and entropy density turns out to be 



■n l a' 

1 = ir + - 

s Air 7i 



c 2 



B 2 



6 +?(9 2 -5r) 



+ 0(a' 2 ). 



(4.18) 



(4.19) 



(4.20) 
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In B — > limit this result matches with [19]. 20 

Let us end up this section with a discussion of the extremal limit. In the absence of 
the moduli, the extremality condition is 2B 2 + rf l (q 2 — 6) = and so the shear viscosity 
to entropy density ratio becomes 



7] 1 a' 
s An n 



B 2 3c 2 / B 



" Cl 7 r+ 2 r r 4 
'h z V 'h 



+ 0(a' 2 ). (4.21) 



Therefore, there is a drastic change when the magnetic field is turned on. That is, 
unlike the electrically charged solution studied in [19], the leading correction of r]/s in 
the extremal limit depends on both, c\ and c 2 . As expected, in B — > limit our result 
matches with the one of [19]. 

In the presence of the scalar field, the extremal limit is more constrained. The 
scalar is fixed by the Chern-Simons coupling only, and by using some of the formulas 
derived in Section 3, namely 

2B 2 e aiph + r A h (q 2 e aiph - 6) = 
q 2 r\ -B 2 = 

^-^-ln|C| =0. (4.22) 
we obtain the following expression for the shear viscosity to entropy density ratio: 

^ = -^(l-8c 1 a , ) + C(«' 2 )- (4.23) 
s An 

Interestingly enough, we observe that there is no dependence of the horizon value of 
the modulus as in the non-extremal case. Also, the coefficient of c 2 vanishes and so the 
dependence of c 2 drops out. 



5. Discussion 

The goal of this paper was two-fold. First, to extend the work of [11] by including 
scalar fields. Second, to compute the low frequency transport coefficients of a (3 + 1)- 
dimensional field theory whose gravity dual is an electrically and magnetically charged 
planar black hole in AdSs Einstein-Maxwell theory with a gauge Chern-Simons term. 

Since, in the extremal case, we have used the entropy function method to study 
the near horizon geometry, we were able to find analytic expressions for the entropy 

20 Note that our 'g' is different than V of [19]. In [19], q is the physical charge (up to some 
normalization). In our case, the physical charge is ~ r\q. 
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density and the near horizon parameters. In the absence of the scalar field we confirm 
the numerical results of [12], in particular the existence of a critical value of the magnetic 
field for which the entropy density vanishes. However, we have found a larger class of 
solutions — we have presented a detailed discussion in Section 3 and we do not want 
to repeat the details here. However, we would like to point out that one of the special 
values of the Chern-Simons coupling, ( = 1, for which there exist finite horizon area 
solutions with both, electric and magnetic, fields non-zero corresponds to a special 
embedding of the gauge field U(l) C SU(4)r [11]. 

In [12], it was proposed that (3 + l)-dimensional gauge theories holographically 
dual to (4 + l)-dimensional Einstein-Maxwell-Chern-Simons theory undergo a quantum 
phase transition in the presence of a finite charge density and magnetic field. The 
authors of [12] also argue that a non-vanishing entropy density is 'exotic' from the 
point of view of CFTs that arise in the AdS/CFT duality. By turning on the magnetic 
field, instabilities can appear and this is consistent also with our analysis. 

However, we believe that this result should be taken with caution. It is well known 
that, in many situations, by adding higher derivative corrections the singularities can 
be 'dressed' by horizons. 21 Also, by turning on the moduli the geometries become more 
stable, though there can appear other kind of instabilities, e.g. the ones associated to 
the AdS2 spacetime. 

Since in string theory, typically, the gauge kinetic terms will couple to various 
scalars [15], we have considered a simple extension of [11] by considering an exponential 
scalar coupling. In the presence of the moduli some important changes occur. First 
of all, in the non-extremal case, the near horizon data depend on the boundary value 
of the scalar. Therefore, the hydrodynamic properties of the dual field theories with 
higher derivative terms will be also controlled by the moduli. 

Let us start with a discussion of the extremal case for which there is a drastic 
change. That is, due to the attractor mechanism, the near horizon geometry is universal 
and does not depend of the asymptotic values of the moduli. The horizon area is finite 
and so there is a large ground state degeneracy. One important question is what is 
happening in the presence of scalars? Is the ground state degeneracy still unstable under 
inclusion of a (large enough) magnetic field? Are there finite horizon area solutions for 
any value of the Chern-Simons coupling when both (magnetic and electric) fields are 
non-zero? 

Before providing concrete answers to these questions, we would like to emphasize 
that the attractor mechanism plays an important role in classifying the bulk theories 

21 In [49], D'Hoker and Kraus have found a new solution with a near horizon geometry that does 
not contain an AdS2- Based on the symmetries, one can argue that higher curvature corrections do 
not change the near horizon form of the solution. 
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where the extremal entropy vanishes. A detailed discussion of the attractor mechanism 
in AdS 5 and its embedding in string theory was presented in [37]. 

The universality of dual theory which means that the IR physics does not depend 
upon the UV details, becomes in the holographic context the statement that the bulk 
solution near the horizon does not depend upon the details of the matter at large values 
of radial coordinate (boundary). Indeed, within the attractor mechanism, the black hole 
horizon (IR region) does not have any memory of the 'initial conditions' (UV values of 
the moduli) at the boundary. This is due to the existence of an infinite throat in the 
extremal near horizon geometry. In the presence of one gauge field (electric or magnetic) 
coupled with one scalar, the extremal limit is a naked (null) singularity. The reason is 
that the effective potential [34] does not have a minimum at the horizon. However, when 
the effective potential has a minimum at the horizon, there exist extremal solutions with 
finite horizon area. In our case we are not able to define an effective potential, but we 
have used, instead, the entropy function formalism. Therefore, the attractor mechanism 
can be regarded as a 'litmus test', which all extremal black holes with finite area should 
pass (the horizon values of the moduli should be finite). 

In the presence of Chern-Simons term, there is no known concrete expression for 
the effective potential (there is a non-trivial fibration for AiS^). Instead, we have used 
the entropy function and shown that the horizon value of the scalar is fixed. What came 
as a surprise was the fact that the modulus is fixed by the Chern-Simons coupling only. 
A simple argument 22 of why is so is that, by rescaling the metric, we can connect 
the (near horizon extremal) solution in the absence of the modulus with the solution 
when the modulus is turned on. The starting point is a finite area solution and so 
we use the solution for £ = ±1 in the absence of the modulus. We then obtain a 
finite area extremal near horizon solution for which the modulus is fixed by the new 
Chern-Simons coupling. Therefore, in the presence of the moduli, there exist finite 
area horizon solutions for any value of the Chern-Simons coupling. We can also safely 
argue that the large ground state degeneracy is stable against scalar field perturbations. 
However, as in the case when the modulus is turned on, there also exists a critical value 
of the magnetic field and so we expect similar physical interpretations as the ones in 
[12]. 

Let us now comment on the second part of our paper. We have computed the 
shear viscosity to entropy density ratio in the presence of the most general bosonic four- 
derivative action (with one electromagnetic field) with moduli dependent couplings. 

In general, there are two types of first-order corrections due to higher derivative 
terms. The entropy/area law is modified due to the additional terms in the action 

22 We would like to thank Rob Myers for a discussion on this point. 
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and/or the modification of the area due to the change of the metric at the horizon (the 
extra terms in the action change the equations of motion). Therefore, one has to use 
Wald formula to compute the entropy. 

To compute the shear viscosity, we have used the method of [17]. Since we needed 
just the near horizon data, the computations are not very complicated even in the 
presence of higher derivative terms. We would like to emphasize that there is no 
ambiguity in defining the overall coefficient of the effective action and the results are 
consistent. Another important observation is that, in the presence of a magnetic field, 
there is a 'preferential' direction in the plasma. A priori it is not clear that the dual 
gravitational mode, h xy , decouples. Indeed, that is not the case when the momentum 
is non-zero — however, in the hydrodynamic limit we have the required decoupling 
[44] and this is what we needed for the computation of the corresponding correlation 
function. 

The shear viscosity to entropy density ratio is controlled by the horizon values of 
the scalars and so, in the non-extremal case, an operator deformation in the QFT will 
produce an interpolating non-trivial flow in which the moduli approach the (IR) black 
hole horizon. In the extremal case, though, the horizon moduli values are fixed and so 
the shear viscosity to entropy density ratio does not depend on the asymptotic values 
of the scalars. Therefore, QFTs with different UV fixed points can flow to the same IR 
fixed point. 

In the presence of the moduli, r]/s also depends only on the c\ (~ (c — a)/c, where 
c and a are the central charges in the CFT) and c 2 (the coupling of the stress tensor 
to the U(l) current) coefficients. The other two coefficients, C3 and C4, parametrize 
couplings of the four-point function of the U(l) currents and so they should appear in 
the expressions of the charge density and conductivity [19]. 

It can be easily checked, that, when the scalars are turned off we recover the results 
of [19]. Since [19] contains a very detailed discussion on the relation with the previous 
results in the literature, we do not want to repeat it here. We would like to point out 
just that, in the extremal limit and in the presence of the magnetic field, we also obtain 
a deconfined 'plasma' in the dual CFT. However, in the presence of the magnetic field, 
rj/s depends on both coefficients, c\ and c 2 — when the magnetic field vanishes, rj/s 
only depends on c 2 and we recover the result of [19]. In N = 2 supergravity, the bulk 
supersymmetry constraints all four derivative couplings to be proportional to a single 
overall constant. 

When the modulus is turned on, the extremal limit is more intriguing. The de- 
pendence of 77/s on c 2 drops out and so rj/s depends only on the central charges in the 
dual CFT. This is somehow similar with the N = 2 supergravity case mentioned above 
for which all the coefficients are fixed by the central charges. 
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We close with some future directions. 

Due to the presence of a Chern-Simons term, there is a more interesting kind of 
instability 23 , which can appear [50]. This instability is caused by a non-normalizable 
mode of the AdS 2 factor in the near horizon geometry. 24 This instability appears for 
large enough Chern-Simons couplings and a finite momentum. 

It will be interesting to check [51] if this kind of instability appears also in our case. 
We would like to point out that the near horizon geometry, when there are background 
electric and magnetic fields, is in fact a fibered AdS2 x R 3 , not a direct product as in 
[50]. In this case an analysis of the near horizon geometry as in our paper (see, also, 
[52]) is more appropriate. 

Using arguments as in [53], a computation of other transport coefficients is also 
possible [14]. However, since the near horizon data are not enough, a more involved 
analysis is needed. Another reason of why the analysis is more complicated is that 
the hydrodynamic limit does not commute with the limit of small magnetic fields. 
Therefore, one expects drastic changes of the transport properties of magnetized fluids. 
In this context, it will be interesting to study how the attenuation and the speed of 
sound waves are affected by the background magnetic field — a similar analysis for a 
(2 + l)-dimensional plasma was presented in [44]. It will also be interesting to see if new 
dissipative transport coefficients can appear in the effective hydrodynamic description 
of plasmas in external magnetic fields [54]. 
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Appendix 

A. Kubo formula and shear viscosity 

In this appendix we present a brief review of Kubo formula in the context of relativistic 
hydrodynamics. Since we are interested in computing the transport coefficients by using 
the AdS/CFT duality, we follow closely the nice reviews [5] (see, also, [55] for a good 
introduction in relativistic hydrodynamics). 

Current understanding is that the matter produced in the heavy-ion experiments 
at RHIC behaves collectively like a fluid. That is, the system rapidly thermalizes and 
comes into /oca/ thermal equilibrium; then it evolves according to hydrodynamics until it 
hadronizes and the temperature becomes smaller than the deconfinement temperature. 

Therefore, the relativistic hydrodynamics — in short, the Navier-Stokes equations 
and their relativistic generalizations — is the best method currently available for mod- 
eling the spacetime evolution of sQGR 

Hydrodynamics can be understood as an effective theory, which describes the dy- 
namics at large length and time-scales. It relies only on the assumption of local equi- 
librium. 25 

In order to describe essential non-equilibrium phenomena in heavy-ion collision, 
a transport theory approach seems necessary. The first important step to investigate 
small perturbations of a high temperature hadronic matter is the linear response for- 
malism. 

If Z m f p is the mean free path — the average distance traversed between collisions 
by particles (of the liquid) — then, a perturbation of the liquid (in equilibrium) will 
disappear at distances of order / m f p . However, the fluctuations associated to some con- 
servation laws can propagate without being damped off. For example, a perturbation 
of (mean) energy can propagate to infinity if there is no internal friction (viscosity) of 
the fluid. 

In the rest of this section we would like to explain why a measure of r]/s is impor- 
tant and we also briefly review the Kubo's formula, which is useful for computing the 
viscosity. 

25 One does not have to make other assumptions on the classical/quantum nature of the phenomena 
involved or the type of particles/fields and their interactions. 
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Unlike in standard thermodynamics where the system is in global thermodynamic 
equilibrium (the intensive parameters, e.g. the pressure (P) and temperature (T), 
are constant throughout the volume) and at rest, we are interested in systems whose 
pressure and temperature vary with space and time, and which are not at rest. 

Though, we request that the system is in local thermodynamic equilibrium. That is, 
the fluid is characterized only by local temperature and velocity fields that vary slowly 
on the scale set by the temperature. Consequently, one can assume the thermodynamic 
equilibrium in some neighborhood about any point. 

In hydrodynamics, it is useful to work with densities per unit volume: the energy 
density e = U/V and the entropy density s = S/V. The usual relation for the energy 
can be rewritten as Ts = e + P — note that all these densities are intensive quantities. 26 

In hydrodynamics, we work directly with the equations of motion and assume an 
expansion in derivatives so that <C l m ] p . We express the stress tensor through the 
temperature T(x) and velocity u^(x). 

For an ideal fluid, the hydrodynamic equations are equivalent with the laws of 
conservation of energy and momentum: 

9^ = 0, T"" = (e + P)u¥ + Pf. (A.l) 

It is clear from the above expression that the momentum density is (e + P)v and so, 
unlike the non-relativistic fluid, the pressure contributes to the inertia of a relativistic 
fluid. 

At the next order we obtain the energy momentum tensor of a relativistic viscous 
fluid whose conservation equation reproduces the relativistic version of the Navier- 
Stokes equation. The stress tensor at this order becomes 

T nv = pg n» + ^ e + P y» u v _ a ^ 

with the dissipative part 



gliv _ p[iapv/3 



(A.2) 



Here, P^ u = g^ v + u^u" is the usual projector operator: for a fluid at rest it becomes 
ptiv _ ( Ji a g(o j +1, +1 5 +1) and so it projects on space. 

The numerical coefficient of the traceless part, 77, is called the shear viscosity. The 
numerical coefficient of the trace part, g, is called bulk viscosity. 



26 All thermodynamic quantities associated with a fluid element, e.g. s, e, and P, arc defined in the 
rest frame and so they are Lorentz scalars by construction. 
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To understand better the physical interpretation of the transport coefficients and 
why a measure of rj/s (and g/s) is relevant (see, also, [56]), let us work in the fluid rest 
frame = (1,0,0,0) and consider a small fluctuation around the thermal equilibrium: 
«* < 1, u° ~ 1 + O(uV), T = T + 5T, etcetera. 

For the ideal fluid, we find the following perturbed stress tensor at first order: 



T uu = e + <5e, T Ul = (e + P)u l , T y = (P + 5P)<f J . (A.3) 

The deviation from the ideal fluid stress tensor (at first order) in the fluid rest 
frame (we use the properties of the projector) is 

STij = (Tij = r] i^djUi + diUj - ^ijdkU 1 ^ + 0Sijd k u k 

= ^ fa? + d t T« - 2 -5 t ,d k T k0 "j + ^d k T k °- (A.4) 

In the last expression we used (A.3) and Ts = e + P. 

Now it is clear that, since the temperature is the only relevant energy scale (ex- 
perimental parameter), the viscous terms are characterized by the coefficients rj/s and 
g/s (medium parameters). The physical interpretation is as follows: the bulk viscos- 
ity encodes the resistance of the system to uniform expansion and the shear viscosity 
controls the rate of the momentum diffusion in the transverse direction to the flow. 

For conformal fluids the energy-momentum tensor should be traceless in flat space, 
which implies g = and e = 3P. 

We have defined the dissipative coefficients as phenomenological constants. How- 
ever, it is important to find a way to compute them directly from the microscopic 
theory. The computation of the shear viscosity in quantum field theory relies on the 
Kubo formula — the dual models of CFTs are well suited to evaluate Kubo expressions. 

Let us now briefly present Kubo formula and its physical interpretation. 

We are interested to understand the response of a macroscopic quantum (many- 
body) system to a localized disturbance, which is created by an applied external force. 
Let us consider a local observable 0(x, t) (e.g., the charge current or T M ", which is a set 
of conserved currents) and an external source that couples linearly to the observable 
so that the new action is 

S = S + J dxO{x)J(x). (A.5) 

The response, which is the change (from zero at the equilibrium) in the expectation 
value of the observable induced by the source, is linear in the source: 

< 0(x) >= J dy X (x - y)J{y) = x ■ J- (A.6) 
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The coefficient of proportionality is called 'susceptibility' and is nothing else than the 
retarded Green function of the physical observable: 

x (x -y) = - l -9{x° - y°) < [0(x), 0{y)} > . (A.7) 
One can also Fourier transform to obtain that, indeed, 

Let us know find the Kubo formula for viscosity (see, e.g., [5]). We consider a 
perturbation of the metric so that the only non-zero component h\2(t) does not depend 
of space coordinates (it corresponds to p = in momentum space). Obviously, this 
perturbation can not excite the temperature (which is a scalar) and the velocity (which 
is a vector). One can easily compute the 12-component of the stress tensor (by using 
the covariant derivative) to obtain 

rj = - lim -Im f / d 4 x6(t)e iult < [T 12 (x), T 12 (0)] > j . (A.9) 
w^O u \J J 

B. Coefficients in expansion (3.3) at 0(r — rh) 2 

For completeness, in this appendix we present the other coefficients that appear in the 
near horizon expansion of different fields in eq. (3.3). They are 



Urf + 3 ) + A ZlTh ~ 

1 / 



v 2 = 



288tt 2 T 2 



^ e -2a<p h ^24B 4 ( 2 e a4,h + 24 (hB 2 + 2q 2 ^j e 3a * h - 72B 2 ( 2 



+ (^6B 4 (a 2 - 4^ - B 2 q 2 (^a 2 + + q 4 ^3a 2 - 4 

W2 = 2 g g ^, 2 ( - 36B 4 C 4 - 245 2 C 2 {b 2 + q 2 ^j e 3a ^ + 48 (q 2 - 2B 2 ^j e^ h + 2885 2 C 2 e 2a ^ 
+ (^-3B 4 (^a 2 - 4j - 8B 2 q 2 + q 4 ^3a 2 - 4^ e 6a<t>h - 144e 4a< ^ 

B( (B 2 (9C 2 e- 3a ^ + l) + 5 (q 2 - 6e- Q ^)) 



6ttT 



Bq(e- 3a * h ((3a 2 + 4) (5 2 - q 2 ) e 3a ^ - 125 2 C 2 + 24e 2a ^) 



Pl 96tt 2 T 2 
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C. Details on entropy function when the modulus is turned on 

In the presence of the modulus, the attractor equation are: 



^V^2 [ AB 2 L 2 e av h + V 2 { _ 4L + 24L 2 + ^2 ^ h + ^2)] = g 

' 7 [4E 2 L 2 e a ^ - w 2 (-4L + 24L 2 + 4gV^ + 3t; 2 z 2 )] = 



2/3tt 



[4B 2 L 2 e a *"' + w 2 (-24L 2 + Aq 2 e aiph + v 2 z 2 )] = 
[L (0 - 45p 2 C) + v 1 ^T 2 {Ap 1 qe aLPh - v 2 zl)] = 
^ [25LC + ^i^ie^ h ] = 



2tt/3 f 85 Pl c — ^-^ — + g ) = 

(OD2 On 2 ' 



(C.l) 



From the last equation, we observe that, in the presence of the modulus, a non- 
trivial extremal black hole can be obtained if neither q nor B vanish. Following the 
same steps as in Section 3.2.2, we can obtain the near horizon geometry. First we 
compute the following quantities: 

e- a ^L(Q + 85piC) 



Zl 



Pi 



4v ly /V2~ 

L{ Pl Q + Q + ABp{Q 



3 



-Q + \/Q 2 ± 16B(2BE^v 2 t ©)C 



8B( 

q 2 - i6ffec 

32e a ^(3v 2 - B 2 e-^(e^ ± 0) 
With these results, the constraint (5th) relation becomes 



v * = — — aVh , 3^ —. (C.2) 



B-^ h /2 ( ^^^.M-B^ 



36V (^±C)A/^ = (C.3) 
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and we obtain the following expression for the entropy: 



y 2e 2 [3v(e 2 — B 2 (e 2 ±Q) 
Now, let us rewrite the near horizon geometry in terms of the charges 



q 6v! V 2 

* = ± , 2B ^ (C.5) 

3 / Q*-KBe £ 



By using the last equation of C.l we obtain 

and from the fifth relation of C.l, we see that the only possibility is a fixed horizon 
value of the scalar: 

¥>fc = §HCI- (C7) 

What comes as a surprise is the fact that the horizon value of the scalar is fixed by the 
Chern-Simons coupling. Consequently, the entropy is 



S = ^\l • ( C - 8 ) 



D. Decoupling of h xy mode 

In what follows, we provide a detailed derivation of the decoupling of the dual gravi- 
tational mode (4.1). We have explicitly checked that the h xy = e lt " +2V ^eQ>(r) mode 
does not couple with any other field when the momentum vanishes. 

For two derivative gravity theory, this can be easily seen from the equations of 
motion (2.3, 2.4, 2.5). However, we are interested in the most general four-derivative 
action (4.9). In this case, instead of writing the equations of motion in the presence of 
higher derivative terms, we will explicitly compute the action up to order e 2 . In this 
way, it can be explicitly checked that there is no coupling between h xy and the other 
fields. 
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Let us turn on the following perturbations of the metric 



(°) , u 



e ift "e£ 2 (r) 
e ltu, e£ A (r) 



e ifc "ef 2 (r) e itu >e£ 3 (r) 
e 2W ^Z(r) 2 - U{r) e ituJ eT(r) 



e itw e£ 4 (r) e^e£ 5 (r) \ 

e 2W(r) Z(r) + e itul ev(r) 



e^eT(r) e 2V » 
e ifcJ+2y ( r )e$(r) 
e 2W/(r) Z(r) + e ittJ et;(r) e ifaJ ex(r) 



,iho+2V(r 



)e$(r) 



e 2V(r) 







e 2W(r) 



and gauge gauge fields 
A a = + ef a 



e Uu >ea r (r), e^ea t {r) - £(r), ^ + e^ea^r), + e itu ea y (r), e itw ea*(r) + P(r) 



(D.2) 



Here and are the background metric (2.8) and the background gauge field; e 
dependent terms are the perturbations. 

Using these field excitations, we can now compute the action. 27 The result is 
complicated, but for our purpose it is enough to pick up the $(r)-dependent terms - 
for concretness, let us write the $(r)-dependent part of the action: 



$(rf 



1 (( c - 2V(r) - w(r) ( U(r)c 2Wlr) (iB 2 c av(r) + AQ(r) 2 e AV{r)+a ^ r) 



-2e w ^U'\r) - Ae 4V ^U'(r) (2V'(r) + W\r)) + e 4V ^ +2W ^ Z\rf + 24e 4V ») 
+2U(r) 2 e 4V(r \2P'(r) 2 e a ^ + e 2W{r \W\r) + AV\r)W\r) + W{rf 

+2W"(r) + 2W\rf - V \r) 2 )) + Uu 2 e AV ^ +2W ^))) 
+$(r)$'(r) (2e 2V( ^ +w ^ (U'(r) + U(r) (3V'(r) + W\r))) 

-{2a'e w ^- 2V ^\-2U{r)V'{r){U{r) (B 2 c 2 + Cl e AV ^ +2W ^ Z'(rf) - 2c 1 U(r) 2 e 4V ^ 
(2V"(r) + W'irf + W\rf) + 2 Cl uj 2 e AV{r) ) + 2c 1 U(r)e mr) U\r) 2 V\r) 

+ Cl e 4y(7 V(r) (2U(r) 2 (V"(r) + 3V'(r) 2 ) + co 2 )))/U(r)) 



27 We use the Mathematica notebook for this computation. We emphasize that our system is sym- 
metric in x~ and y— directions. 
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+$(r)$"( r ) 



(2U(r)e 2V ^ +w ^ - 4c ia 'U(r)e 2V ^ +w ^(U'(r)V'(r) 



+2U(r) (V"(r) + V'(r) 2 ))) 



+$ (r 



'^2 



ck e 



/_W(r)-2V(r) 



(f/(r)( J B 2 c 2 + 2c 1 e 4y Ml7'(r)\/ , (r) - c x e w ^+ 2W ^Z\r)< 



+ Cl e 4V » {U\rf + Au 2 ) + 2c 1 f/(r) 2 e 4V ' (r) (-2F"(r) + V\r) 2 + W(r) 2 )) 

+-[/(r)e 2F W +Mf W) 
2 



+$'(r)$"(r) 



(2c 1 a'f/(r)e 2V ' (r)+PF(r) (C/'(r) + 4f/(r)V'(r))) 



+$"(r) 2 



2 Cl a'f/(r) 2 e 2y M +H/ W 



(D.3) 



Since there is no mixing between $ and the other modes in the zero momentum 
limit, this mode remains massless and a computation of the related Green's function 
from the near-horizon data is still possible. 

However, the other metric fluctuations namely h xz or h yz are not decoupled from 
the gauge field perturbations in the zero momentum limit (even in the absence of higher 
derivative terms) — these coupled terms are proportionl to B or ( (CS term). There- 
fore, these modes are not massless in the zero momentum limit and the near horizon 
geometry is not sufficient to compute the related two point correlation functions, i.e. 
< [T XZ ,T XZ ] > or < [T yz ,T yz ] > . 
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